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Abstract 

We prove that the sequence (l/F„+2)n>o of reciprocals of the Fibonacci numbers is 
a moment sequence of a certain discrete probabihty, and we identify the orthogonal 
polynomials as little g-Jacobi polynomials with q = {1 — ^/5)/{l + \/5). We prove 
that the corresponding kernel polynomials have integer coefficients, and from this 
we deduce that the inverse of the corresponding Hankel matrices have 
integer entries. We prove analogous results for the Hilbert matrices. 
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1 Introduction 

In [12] Richardson noticed that the Filbert matrices 

= , 0<t,j<n, n = 0,l,..., (1) 

where F„, n > is the sequence of Fibonacci numbers, have the property that 
all elements of the inverse matrices are integers. The corresponding property 
for the Hilbert matrices {l/{i + j + 1)) has been known for a long time, see 
Choi [5]. Richardson gave an explicit formula for the elements of the inverse 
matrices and proved it using computer algebra. The formula shows a remark- 
able analogy with Choi's corresponding formula for the elements of the inverse 
Hilbert matrices in the sense that one shall replace some binomial coefficients 
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by the analogous Fibonomial coefficients 



with the usual convention that empty products are defined as 1. These co- 
efficients are defined and studied in [10] and are integers. The sequence of 
Fibonacci numbers is Fq = 0, Fi = 1, . . . , with the recursion formula Fn+i — 
Fn + Fn-i, n > 1. 

The Hilbert matrices are the Hankel matrices (sj+j) corresponding to the 
moment sequence 

Jo 

and that the reciprocal matrices have integer entries can easily be explained 
by the fact the corresponding orthogonal polynomials, namely the Legendre 
polynomials, have integer coefficients. See section 4 for details. 

The purpose of the present paper is to show that (l/F„+2)n>o is the moment 
sequence of a certain discrete probability. Although this is a simple conse- 
quence of Binet's formula for F„, it does not seem to have been noticed in the 
literature, cf. [11]. We find the corresponding probability measure to be 

oo 

/^ = (1-?')E?'''^.V^' (3) 

fe=0 

where we use the notation 

1+V5 1- 75 1 

and 6a denotes the probability measure with mass 1 at the point a. The number 
is called the golden ratio. 

The corresponding orthogonal polynomials are little g-Jacobi polynomials 

Pn{x;a,b;q) = 20i | ■,q,xq\ , (5) 

aq 

see [6], speciaUzed to the parameters a = q,b = 1, with q taking the value 
from (4). 

To be precise we define 

Pn(x) := F„+ip„(x0; q, 1; q), (6) 



2 



and these polynomials have integer coefficients, since they can be written 

n 
A;=0 

The orthonormal polynomials with respect to // and having positive leading 
coefficients are given as 

Pn{x) = (-l)("^')^F2„+2Pn(a;), (8) 
so the kernel polynomial 

n 

K^{x,y)^^Pk{x)Pk{y), 

k=0 

is a polynomial in x, y with integer coefficients. If we denote a^^ the coefficient 
to x''y^ in the kernel polynomial, then it is a general fact that the matrix 

^n = (aS^ 0<i,j<n (9) 

is the inverse of the Hankel matrix of the problem (sj+j)Q, see Theorem 2.1 
below. 

This explains that the elements of the inverse of the matrix (Fj_|_j_|_2)o are inte- 
gers, and we derive a formula for the entries from the orthogonal polynomials. 

The Filbert matrices (1) are not positive definite but non-singular, and they 
are the Hankel matrices of the moments of a (real- valued) signed measure with 
total mass 1. The orthogonal polynomials for this signed measure are the little 
g-Jacobi polynomials 

Pnix<P;l,l;q) = E(-l)'^"-© Q.Cn')/^ (10) 

fc=0 

and a simple modification of the positive definite case leads to Richardson's 
formula for the entries of the inverse of the Filbert matrices. 

The two results can be unified in the statement that for each o; e N = 
{1,2, . . .} the sequence (F„/Fa+,„)„>o is a moment sequence of a real-valued 
measure /la with total mass 1. It is a positive measure when a is even, but a 
signed measure when a is odd. The orthogonal polynomials are little g-Jacobi 
polynomials Pnix(f); q°'~^, 1; q)- This is proved in section 3. 

In section 2 we recall some basic things about orthogonal polynomials both 
in the positive definite and in the quasi-definite case, and Theorem 2.1 about 
the inverse of the Hankel matrices is proved. 
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In section 4 we briefly discuss the matrices (l/(« + i + j))o, where a > 0. 
They are related to Jacobi polynomials transfered to the interval ]0,1[ and 
belonging to the parameters (0, a — 1). This leads to a generalization of Choi's 
result, which corresponds to a = 1. 

After the circulation of a preliminary version of this paper Ismail has extended 
the results of section 3 to a one parameter generalization of the Fibonacci 
numbers, cf. [9]. 



2 Orthogonal Polynomials 



We start by recalling some simple facts from the theory of orthogonal polyno- 
mials, cf. [1] or [8] and in particular [4] for the quasi-definite case. 



The positive definite case. 

We consider the set A^* of probability measures on R with moments of any 
order and with infinite support. The moment sequence of e A^* is 



Sn = Snil^) = y" ^" dn{x), n = 0, 1, . . . , 



and the corresponding Hankel matrices are given by 



^So Si 
Sl S2 



S ^ 

S2n J 



, n = 0, 1,... 



:i2) 



The orthonormal polynomials (P„) for /i are uniquely determined by the equa- 
tions 

j Pn{x)Pm{x) dll{x) = 6n,m: n, m > 0, (13) 

and the requirement that P„ is a polynomial of degree n with positive leading 
coefficient. This coefficient is equal to 



Dn-l/ Dni 



(14) 



where D„ = det The reproducing kernel for the polynomials of degree < n 
is defined as 



K^{x,y) = Y.Pu{x)Pk{y), 

k=0 



(15) 
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and is called the kernel polynomial. It is clear that we can write 



n n 



i^n(x,y) = EE«SW. (16) 

i=0 j=0 

where the numbers a^^j are uniquely determined and satisfy afj — a^^^ . If we 

collect these numbers in an (n + 1) x (n + l)-matrix A„ = {^t^j)^ then it is the 
inverse of the Hankel matrix 



Theorem 2.1 



where En is the unit matrix of order n+ 1. 
Proof. For < k < n we have 

k 

/ x''Kn{x,y)dfj.{x) = J2 Prn{y) I )d^i{x), (17) 

which is a polynomial in y of degree k. On the other hand we have 

/n n 
x'^Knix, y) dii{x) = E(E 
j=0 i=0 

and therefore 

n 

when k < j < n, and when j = k the sum equals the coefficient to y'^ in (17), 
i.e. equal to 

\/Dk-i/Dk I x^Pkix) dii{x) = j Pl{x) dii{x) = 1. 

Since the matrix H^An is symmetric, the above shows that it equals the unit 
matrix. □ 



The quasi- definite case. 

If /i is a real-valued signed measure on M with total mass 1 and moments 
of any order, one can still define the moments (11) and the corresponding 
Hankel matrices (12). To define orthogonal polynomials one has to assume 
that (12) is a non-singular matrix for any n, i.e. that the determinants sat- 
isfy Dn = det Hn 7^ 0. On the other hand, if orthogonal polynomials exist 
with respect to a signed measure, then the Hankel determinants are non- 
zero. See [4, Theorem 3.1] for details. In this case the orthonormal polyno- 
mial Pn is uniquely determined by the requirement that the leading coefficient 
J Dn-i/Dn is either positive or purely imaginary with positive imaginary part. 
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The corresponding kernel polynomial Kn has real coefficients, and Theorem 
2.1 remains valid. 



3 Fibonacci numbers 



The Fibonacci numbers can be given by the formula 

F„ = i=(</.--0"), n>Q (18) 

usually called Binet's formula, but it is actually older, see [10], [11]. Here 
, 1 + ^5 ^ 1-v^ ^ , 



Using the number q — (f)/4>, satisfying — 1 < g < and already defined in (4), 
leads to 

F„ = i=0"(l-g"), #^ = -1, (19) 
and for a G N and n > 



oo 



^ a+n V ^ H k=o 

which is the n'th moment of the real-valued measure 

oo 

/^a = (l-?")E?"''^9V0 (20) 

k=0 

with total mass 1. When a is even then /Xq, is a probability measure, but when 
a is odd the masses q°''^ change sign with the parity of k. Note that /i2 is the 
measure considered in (3). 

For the Fibonomial coefficients defined in (2) one has 
and they satisfy a recursion formula 

= ^^-^tA + ^ > ^ > 1, (21) 

see [10], which shows that the Fibonomial coefficients are integers. From (2) 
it is also clear that 

~ {n-k}p-> < k < n. 
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In [6, Section 7.3] one finds a discussion of the little g-Jacobi polynomials 
defined in (5), and it is proved that 

f:Pn{q';a,b;q)Pm{q';a,b;q)^^^{aq)' = (22) 
^ [q;q)k hn[a,b;q) 

where 

^ / ^. x {abq; q)n{l - abq^''+^){aq; q)n{aq; q)oo , .^s) 
{q]q)n{^ - abq){bq]q)n{abq^]q)oo 

In [6] it is assumed that < q,aq < 1, but the derivation shows that it 
holds for |g| < 1, \a\ < 1, \b\ < 1, in particular in the case of interest here: 
— l<q<0,a = q"'^^, 6 = 1, in the case of which we get 

EPn(g^ q--\ 1; q)Pm{q'; q^-', 1; q)q-' = ^n,^ _%n) ' (24) 

This shows that the polynomials 

p„(a;0; 1; g) 
are orthogonal with respect to and that 

j Pnixcj); 1; q)Pm{x(t)] 1; q) dii^{x) = ^n,m ^^^^ _ ^^a+ln^ ■ 

To simplify this apply (19) to get 
(l-?")g°"(g;g) 

n / -\\an Fg I TT Fi+j \ / -\\an Fg fa+n—l\ 



\j=0 



Theorem 3.1 Let a e N. The polynomials p^^^x) defined by 

P^:\x) = {"^:-')^^{x<P;q'^-\l;q) (25) 

can be written 

n 

p^:\^) = IK-^"'^ (26) 

fe=0 

and they satisfy 

j P^:\^)pt\^) df^aix) = 5„,^(-l)-^, (27) 
so the corresponding orthonormal polynomials are 

P^-^ (X) = ^(-l)-F„+2n/i^aPi"^ (X) ■ (28) 
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Proof. By definition, see (5) 



k=o 



/a+n-l\ 



fe=0 



{q(t>xy 



where 



is the g'- binomial coefficient. Using (19) leads to 



n ' 



hence 



k=0 



fc-1 

n 

j=0 



F I • 



which by (19) can be reduced to (26). □ 



Remcirk 3.2 The polynomials p^n\'^) for a = 1 and a — 2 are the polyno- 
mials in (10) and in (7) respectively. 

Corollary 3.3 For a e N 

det(l/F„+,+,)^= ((-l)"mF.nF.^,,("+2^'=-i)^^ \ 



k=l 



which is the reciprocal of an integer. 

Proof. From the general theory it is known that the leading coefficient of the 



orthonormal polynomial P^^^ is ^ Dn-i/ D^, where 
From (26) and (28) we then get 

hence 

J^k k=i 

and the formula follows. □ 
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Theorem 3.4 The i,j'th entry of the inverse of the matrix {l/Fa+i+j)Q is 
given as 

V / o'+^+jy n-j Jj^y n-i Jf\ i Jj^y j \ > 



Proof. From Theorem 2.1 we get 

where a^l^j{a) is the coefficient to x'^y^ in the kernel polynomial Kn{x,y) for 
the orthonormal polynomials -P^"-*. Inserting the expressions (26) and (28) in 
the kernel polynomial and changing the order of summation gives 

n 

k=ui3x{i,j) 

where we for A; > i, j have defined 

C^^Hk;i,j) {-l)''''^''-^'-^^-^^F^^,,(;i^^^^^ 

(30) 

To prove that this expression can be summed to give (29), we use induction in 
n. By symmetry we can always assume i > j. The starting step n = k = i > j 
is easy and is left to the reader. For the induction step let R^°'\n] i,j) denote 
the expression (29). It has to be established that 

(n + 1; i, j) - (n; i, j) = C^") (n + 1; i, j) . 

The left-hand side of this expression can be written 

(_l\in+l){a+i+j)-[l)-{>^)p . (a+i+j-l\ fa+i+j-l\ j, 

where 

V n+l-j J-f \ n+l-i Jw ^ ^ \ n-j /pV n-i Jf 

{Fg+n+i ' ' ' -Pg+i+j+l) (-Pg+n+j ' ' ' -^g+i+j+l) 

(-^1 • • • Fn+-i-j){Fi ■ ■ • Fn+i-i) 

Fa+n+i+lFa+n+j+1 ~ {~^)"~^^'^'' Fn+l-iFn+i-j ■ 

By Lemma 3.5 below (with n replaced by n + 1), the expression in brackets 
equals Fa+2n+2Fa+i+j , and now it is easy to complete the proof. □ 

Lemma 3.5 For n > i, j > and a > the following formula holds 

Fa+2nFa+i+j -^g+n+i-^a+n+j? ( 1) ''-^n— i-^n— ("^-^) 
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Proof. Using Binet's formula, the right-hand side of (31) multiphed with 5 
equals 

Using (fxf) — —1 one gets after some simplification 

which establishes the formula. □ 

Remark 3.6 For a — 1 the expression (29) reduces to 

^_l)n(i+J+l)-(^)-(i)p /n+i+l\ (n+j+l\ /i+jV 

which is the expression found by Richardson [12], except that he expressed 
the sign in a different but equivalent manner. 



4 The Hilbert matrices 

For a > the matrices 

= («/(« + i+i))^, n = 0,l,..., (32) 
are the Hankel matrices for the moment sequence 

s^^) = a t x'^x'^-^dx^-^, n = 0,l,... 
Jo a + n 

of the measure cTq = aa;"~^l]o,i[(a;) dx. The corresponding orthogonal polyno- 
mials are easily seen to be 

ri»){x) = -x-'^+'D- [x^-'+-{l - xr] = (-1)" X: (l) ("Tl - 1)'^""', 
^- k=0 

(33) 

since they are Jacobi polynomials transfered to ]0, 1[, cf. [2]. Using the binomial 
formula for (x — 1)*^ we find 

n 
3=0 

where 

10 



%=Ee)G)("T")=i:ej%';,)("^: 

k=j 1=0 



fn\fa-l+n\ jp ( -n + j, -n - a + j + 1 \ / n\ f2n+a-j-l\ 
^Ul j )^^\ j + 1 '^J = Wl n j' 

where the is summed by the Chu-Vandermonde formula, cf. [2, p. 67]. This 
gives 

n 

ri"H^) = E(-l)^G)("^t''">^- (34) 

j=0 

The orthonormal polynomials with positive leading coefficients are given as 



V a 

so the corresponding kernel polynomials have coefficients a\^j{a) which by 
Theorem 2.1 satisfy 

= ± (« + 2fc)©e)('^+\+-^)("+\+-^). (35) 

k=max 

Theorem 4.1 The i, j 'th element of the inverse matrix of {l/{a + i + j))q is 
given as 

{-iy+^{a ("+:^+*) ("+::+^) ("+^+^-^) . (36) 

In particular they are integers for a Furthermore, 

det{l/{a + i+j)r,^ (^af[{a + 2k){''^f-'y^ . (37) 

Proof. Let R{n;i,j) denote the number given in (36), and define 

C{k;i,j) = i-iy^^ia + 2fc)g) g) ("+\+^-^), k > ij. 

We shall prove that 

n 

R{n;i,j)= J2 C{k;i,j) 

fc=max(jj) 

by induction in n and can assume i > j. This is easy for n = A; = i and we 
shall establish 

R{n + l;i,j) - R{n;i,j) = C{n + l;i,j). (38) 
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The left-hand side of this expression can be written 
where 

_ / a+n+l+i\ / a+n+l+j\ _ / a+n+i\ / a+n+j\ 
\ n+l-j J \ n+l-i J \ n-j J \ n-i J 

^ {{a+n+i)-{a+i+j+l)){{a+n+j)-{a+i+j+l)) _ 
(n+l-j)!(n+l-i)! 

[{a + n + l + i){a + n + 1 + j) - (n + 1 - j){n + 1 - i)]. 

The quantity in brackets equals (a + 2n + 2) (a + i + j), and now it is easy to 
complete the proof of (38). 

The leading coefficient of R!;°'\x) is 



Dn-i _ a + 2n a + 2n-l\ 

where 

Dn = det {a/ {a + i + j))^ = q;"+^ det {l/{a + i + j))^ . 

Therefore 

-^n k=l « ^^-^ 

which proves (37). □ 

Replacing a; by 1 — a;, we see that r^"-*(l — x) are orthogonal polynomials with 
respect to the probability measure a{l — a;)"~^l]o,i[(a;) dx. The corresponding 
moment sequence is 

= -p^' (39) 



and the corresponding orthonormal polynomials are ^J{(y.-\- 2n)/a r^°'^(l — x). 
Therefore 

K^{x, y)^± ^^rt\l - x)rt\l - y), (40) 



k=0 ^ 



showing that the coefficient to x^y^ in aKn{x,y) is an integer when a e N. 
This yields 

Theorem 4.2 Let a eN. The inverse of the matrix 

(41) 



a 







has integer entries. 
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It is not difficult to prove that 



ri"Hi-^) = E(-ir"'©(' 



a+n+fc— 1 
k 



) 




fc=0 



and it follows that the entries of the inverse of (41) are given as 



n 



)■ 



k=max{i,j) 



This formula holds of course for any a > 0. 

The results of this section for a = 1,2 have been treated in the survey paper 
[3] , written in Danish. For a = 1 the formula for the elements of the inverse of 
'H^"-' was given in [5], while the formula for its determinant goes back to Hilbert 
in [7] . In this case the polynomials r^^^ (x) are the Legendre polynomials for the 
interval [0, 1], cf. [2, Section 7.7]. These polynomials have succesfuUy been used 
in the proof of the irrationality of C(3) . For a = 2 we have (a + 2k) /a = 1 + k, 
so the coefficient to x'^y^ in (40) is an integer. In this case Theorem 4.2 can 
be sharpened: The inverse of the matrix (i/(^^2^"'))q has integer coefficients. 
This result is also given in [12]. 
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Abstract 

We prove that the sequence (l/F„+2)n>o of reciprocals of the Fi- 
bonacci numbers is a moment sequence of a certain discrete probabil- 
ity, and we identify the orthogonal polynomials as little g-Jacobi poly- 
nomials with q = {1 — \/5) /(! + \/5). We prove that the corresponding 
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that the inverse of the corresponding Hankel matrices (l/Fj+j_|_2) have 
integer entries. We prove analogous results for the Hilbert matrices. 
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1 Introduction 

In [H] Richardson noticed that the Filbert matrices 

= , 0<z,j<n, n = 0,l,..., (1) 

where n > is the sequence of Fibonacci numbers, have the property that 
all elements of the inverse matrices are integers. The corresponding property 
for the Hilbert matrices {l/{i + j + 1)) has been known for a long time, see 
[Z] 5 [6] , [15] . The last reference contains a table of the inverse Hilbert matrices 
up to n = 9. 

Richardson gave an explicit formula for the elements of the inverse Fil- 
bert matrices and proved it using computer algebra. The formula shows a 
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remarkable analogy with the corresponding formula for the elements of the 
inverse Hilbert matrices in the sense that one shall replace some binomial 
coefficients (^) by the analogous Fihonomial coefficients 

with the usual convention that empty products are defined as 1. These 
coefficients are defined and studied in [12j and are integers. The sequence 
of Fibonacci numbers is Fq = 0,Fi = 1,..., with the recursion formula 
Fn+i = F^ + Fn^i, n>l. 

The Hilbert matrices are the Hankel matrices (sj+j ) corresponding to the 
moment sequence 

Sn = l/(n + 1) = / x" dx, 



JO 

and that the reciprocal matrices have integer entries can easily be explained 
by the fact the corresponding orthogonal polynomials, namely the Legendre 
polynomials, have integer coefficients. See section 4 for details. 

The purpose of the present paper is to show that (l/F„+2)n>o is the 
moment sequence of a certain discrete probability. Although this is a simple 
consequence of Binet's formula for F„, it does not seem to have been noticed 
in the literature, cf. [13]. We find the corresponding probability measure to 
be 

oo 

where we use the notation 

and Sa denotes the probability measure with mass 1 at the point a. The 
number is called the golden ratio. 

The corresponding orthogonal polynomials are little g- Jacobi polynomials 

Pn(x;a,6;g) = 201 I ;q,xq\, (5) 

see [8], specialized to the parameters a = q,b = 1, with q taking the value 
from (Hi). 
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To be precise we define 

Pnix) := F„+ip„(x0; q, 1; q), (6) 
and these polynomials have integer coefficients, since they can be written 

n 

Pn{x)=Y^{-ir-(^{xrt^\x^. (7) 

A;=0 

The orthonormal polynomials with respect to /i and having positive leading 
coefficients are given as 

Pnix) = (-1)("^')/FW^p„(x), (8) 
so the kernel polynomial 

n 

K^ix,y) = Y,Pk{x)Pkiy), 

is a polynomial in x, y with integer coefficients. If we denote a-""* the co- 
efficient to x'^y^ in the kernel polynomial, then it is a general fact that the 
matrix 

^n = (ag), 0<t,3<n (9) 

is the inverse of the Hankel matrix of the problem {si^j)Q, see Theorem 12.11 
below. 

This explains that the elements of the inverse of the matrix (1/^44.^+2)0 
are integers, and we derive a formula for the entries from the orthogonal 
polynomials. 

The Filbert matrices ([1]) are not positive definite but non-singular, and 
they are the Hankel matrices of the moments of a (real-valued) signed mea- 
sure with total mass 1. The orthogonal polynomials for this signed measure 
are the little g-Jacobi polynomials 

n 

p„(x0; 1, 1; q) = {1),^%^'^ (10) 

fc=0 

and a simple modification of the positive definite case leads to Richardson's 
formula for the entries of the inverse of the Filbert matrices. 
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The two results can be unified in the statement that for each a G N = 
{1, 2, . . .} the sequence {Fa/ Fa^n)n>Q is a moment sequence of a real- valued 
measure with total mass 1. It is a positive measure when a is even, 
but a signed measure when a is odd. The orthogonal polynomials are little 
g-Jacobi polynomials Pn{x4>; q^~^ , 1; q)- This is proved in section 3. 

In section 2 we recall some basic things about orthogonal polynomials 
both in the positive definite and in the quasi-definite case, and Theorem 12.11 
about the inverse of the Hankel matrices is proved. 

In section 4 we briefly discuss the matrices (l/(a + ^ + j))o, where a > 0. 
They are related to Jacobi polynomials transfered to the interval ]0, 1[ and 
belonging to the parameters (0, a — 1). This leads to the formula fl36p . which 
for a = 1 is the formula for the elements of the inverse Hilbert matrices. 

After the circulation of a preliminary version of this paper (dated April 
10, 2006), Ismail has extended the results of section 3 to a one parameter 
generalization of the Fibonacci numbers, cf. |11] . 



2 Orthogonal Polynomials 



We start by recalling some simple facts from the theory of orthogonal poly- 
nomials, cf. [1] or [To] and in particular [5] for the quasi-definite case. 

The positive definite case. 

We consider the set A^* of probability measures on R with moments of 
any order and with infinite support. The moment sequence of G A^* is 



Sn = Sn(Ai) = / dfj,{x), n = 0, 1, 



and the corresponding Hankel matrices are given by 



( So Si 



Sn \ 

Sn+1 

S2n ) 



, ra = 0, 1, 



(12) 



The orthonormal polynomials (P„) for yU are uniquely determined by the 
equations 

P„(x)P„(a;) d[i{x) = 6n,m, n,m>0, (13) 
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and the requirement that P„ is a polynomial of degree n with positive leading 
coefficient. This coefficient is equal to 

^jDn-l/Dn, (14) 

where D.^ = dctif„. The reproducing kernel for the polynomials of degree 
< n is defined as 

n 

J]Pfe(x)Pfe(y), (15) 

fc=0 

and is called the kernel polynomial. It is clear that we can write 

11 n 
i=0 j=0 

where the numbers afj are uniquely determined and satisfy a^/j = a^J^} . If 

we collect these numbers in an (n + 1) x (n + l)-matrix = (^1 j'*), then it 
is the inverse of the Hankel matrix 



Theorem 2.1 



A H — H A — T 



where /„ is the unit matrix of order n + 1. 

Proof. For < k < n we have by the reproducing property 

J x''K„{x,y)dfi{x)^y\ (17) 

On the other hand we have 

n n 

x^K^{x,y) dn{x) = J^(J^Sfe+iaJ^.^)y^ 

j=0 1=0 

and therefore 

n 

(i) r 
2_^^k+i%j - dk,j- 

1=0 

□ 

The quasi- definite case. 
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If is a real-valued signed measure on M with total mass 1 and moments 
of any order, one can still define the moments fllip and the corresponding 
Hankel matrices f lT2|) . To define orthogonal polynomials one has to assume 
that f|T2|) is a non-singular matrix for any n, i.e. that the determinants satisfy 
Dn = det Hn 7^ 0. On the other hand, if orthogonal polynomials exist with 
respect to a signed measure, then the Hankel determinants are non-zero. 
See P, Theorem 3.1] for details. In this case the ortho normal polynomial 
Pn is uniquely determined by the requirement that the leading coefficient 
\/ Dn-\l Dn is either positive or purely imaginary with positive imaginary 
part. The corresponding kernel polynomial has real coefficients, and 
Theorem 12.11 remains valid. 

3 Fibonacci numbers 

The Fibonacci numbers can be given by the formula 

F„ = i=(0"-0"), n>0 (18) 



usually called Binet's formula, but it is actually older, see [12], [IS]- Here 

, l + VS ^ 1-V5 



2 ' ^ 2 



1 - 



Using the number q = satisfying — 1 < g < and already defined in 
(HD, leads to 

F„ = -^0"(l-g"), q<f>' = -l, (19) 



and for a; G N and n > 

Fa V5Fa 1 



1 - 



which is the n'th moment of the real-valued measure 

oo 

fia = {l-qnJ2^''%'/<^ (20) 



fe=0 
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with total mass 1. When a is even then probabihty measure, but 

when a is odd the masses g"'' change sign with the parity of k. Note that fi2 
is the measure considered in ([3]). 

For the Fibonomial coefficients defined in ([2]) one has 

□ f = 1, < < n < 2, 

and they satisfy a recursion formula 

a)F=^^-i(V)F+^-^+i(r;)r ^>^>i> (21) 

see [12], which shows that the Fibonomial coefficients are integers. From ([2]) 
it is also clear that 

(l)^ < A; < n. 

Vfc/F Vn— fc/F' — — 

In [HI Section 7.3] one finds a discussion of the little g-Jacobi polynomials 
defined in ([5]), and it is proved that 

£p„(g^ a, 6; q)pUq'; a, 5; g)^^(ag)'= = , f"''? (22) 

where 

/i (a = (^^g- g)^(l ~ g)n(Qg; g)oo /^^\~« (-23) 

{q;q)n{l - abq){bq;q)n{abq'^;q)oo 

In [8] it is assumed that < q,aq < 1, but the derivation shows that it 
holds for |g| < l,|a| < < 1, in particular in the case of interest here: 
— 1 < g < 0, a = g°~^, 6 = 1, in the case of which we get 

J2p.{q'; q--\ 1; g)p„(/; 1; q)q-' = ^n.^ (i _ g"^^^.) • (24) 

This shows that the polynomials 

are orthogonal with respect to and that 

.c-i . . . (1 - g")g°"(g; g)^ 



g" , 1; q)pm{x(j); g" , 1; g) = 5„,^ 



(g-g)2(l-g-+2n)- 



To simplify this apply (IT^ to get 



\j=0 



Theorem 3.1 Let a E N. The polynomials pn (x) defined by 
p(^\x) = {''+:^')^Pn{xc^-q--\l-q) 

can be written 



k=0 

and they satisfy 

so the corresponding orthonormal polynomials are 

PjC\x) = ^J{-lY^F^+2n/FJ:\x). 

Proof. By definition, see ([5]) 

fe=0 ^^'^ 



fc=0 

where 



1 {<i;<l)k{q;q)n-k 
is the g-binomial coefficient. Using (1T9|) leads to 



n 






A; 


9 






F 



hence 



k-l 



fc=0 j=0 

which by ( |T9i) can be reduced to f l26l) . □ 
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Remark 3.2 The polynomials pn \x) for a = 1 and a = 2 are the polyno- 
mials in ffTOj) and in ([7]) respectively. 

Corollary 3.3 For a e N 

det(i/F„+,+,)^ = F«nF«+2.rr')F 

V k=l 

which is the reciprocal of an integer. 

Proof. From the general theory it is known that the leading coefficient of 
the orthonormal polynomial Pi°^ is ^/1\~\J1\^ where 

Dn = det(Fa/Fa+j+j)Q. 

From (!26l) and (p8l) we then get 
hence 

" fc=l k=l 
and the formula follows. □ 

Theorem 3.4 T/ie i, j 't/i entrj/ o/ ^/ie inverse of the matrix {1 / Fa+i+j)^ is 
given as 

{-l)<'^^^-^^H^<i) F^^^^.i^:^^^^^ (29) 
Proof. From Theorem 12.11 we get 

where a^l^j{a) is the coefficient to x^y^ in the kernel polynomial Kn{x,y) for 

the orthonormal polynomials Pn"\ Inserting the expressions fl2^ and f l25]) 
in the kernel polynomial and changing the order of summation gives 

n 

fc=max(i j) 

9 




where we for k > i,j have defined 

(30) 

To prove that this expression can be summed to give we use induction 
in n. By symmetry we can always assume i > j. The starting step n = k = 
i > j is easy and is left to the reader. For the induction step let R^"\n; 
denote the expression fl29|) . It has to be established that 

(n + 1; j) - {n-z,j) = C^") (n + 1; 2, j) • 
The left-hand side of this expression can be written 

where 

IP /a+n+l+A /a+n+l+j\ / /a+n+iN /a+n+j\ 

\ n+l—j /jj-V n+l—i /W ^ ^ V n~j Jf\ n—i Jf 

_ i^a+n+i ■ ■ ■ -^a+j+j+l)(-^a+r!.+j ' ' ' -^a+i+j+l) 
(Fi ■ ■ ■ Fn+l-j){Fi ■ ■ ■ Fn+l-i) 

[Fa+n+i+lFa+n+j+1 ~ (~l)°^*^''-^n+l-j-^n+l-j] ■ 

By Lemma [3.51 below (with n replaced by n + 1), the expression in brackets 
equals Fa+2n+2Fa+i+j, and now it is easy to complete the proof. □ 

Lemma 3.5 For n > i, j > and a > the following formula holds 

Fa+2nFa-\'i-\-j -Fa+n+i-^a+n+j ( 1) Fn~iFn~-j ■ (31) 

Proof. Using Binet's formula, the right-hand side of (13T|) multiplied with 
5 equals 

Using 00 = —1 one gets after some simplification 

which establishes the formula. □ 

Remark 3.6 For a = 1 the expression ( l29l) reduces to 

V ^) -'J+J+IV n-j n-i i Iw' 

which is the expression found by Richardson [H], except that he expressed 
the sign in a different but equivalent manner. 
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4 The Hilbert matrices 



For a > the matrices 

n^^) = {a/{a + i + 3))l, n = 0,l,..., (32) 
are the Hankel matrices for the moment sequence 

= a / dx = n = 0, 1, . . . 

Jo a + n 

of the measure cTq = ax°~^l]o,i[(a;) rfx. The corresponding orthogonal poly- 
nomials are easily seen to be 



1 " 



(33) 

since they are Jacobi polynomials transferee! to ]0, 1[, cf. |3]. Using the 
binomial formula for (x — 1)'^ we find 



-l)-5^(-l)^x"-^c,, 

j=0 



where 



9 = EC)a)m=E(7)(;j("7ir) 

/'raWa-l+ra\ 17 /-n + i, -n - a + j + 1. i ^ {n\ f2n+a-j-l\ 

where the 2-^1 is summed by the Chu-Vandermonde formula, cf. [3, p. 67]. 
This gives 

n 
j=0 

The orthonormal polynomials with positive leading coefficients are given as 
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so the corresponding kernel polynomials have coefficients a!i'^{Q) which by 
Theorem 12.11 satisfy 



= (-1)^-^^- (a + 2A:)(5) g) ("+\+^-^) ^'t^-'). (35) 

fc=max (i j) 

Theorem 4.1 The i,j 'th element of the inverse matrix o/ (l/(« + i + j))q 
is given as 

{-lY+^{a + ("n"?) ("^?"') ("^^f ) • (36) 

In particular they are integers for a G N. Furthermore, 

det(l/(a + z + j))o= i^^{a + 2k){^+'^-'f^ . (37) 

Proof. Let i?(n; z, j) denote the number given in fl5^ . and define 
C(A:; ^, j) = (-l)^+^ (a + 2k) (1) r'^^^'^^), k > 

We shall prove that 

n 

R{n;i,j)= '^(^'^'^') 

fc=max{i j) 

by induction in n and can assume i > j. This is easy for n = k = i and we 
shall establish 

i?(n + 1; 2, j) - R{n; t,j) = C{n + l;t,j). (38) 
The left-hand side of this expression can be written 

where 



V n+l—j J V n+1— j / V n—j ) \ n—i ) 

^ {{a+n+i)---{a+i+j+l)){(a+n+j)-{a+i+j+l)) _ 
(n+l-j)!(n+l-i)! 

[{a + n + l + i){a + n + l+ i) -{n + l- i){n + l- i)]. 
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The quantity in brackets equals {a + 2n + 2){a + i + j), and now it is easy 
to complete the proof of 

The leading coefficient of is 



Dn-i a + 2n (a + 2n — 1 



Dn V a V n 



where 

Dn = det {a/ (a + i + j))^ = det + i + j))^ . 

Therefore 



fc=i fe=i 



which proves (l37j) . □ 

Replacing x by 1 — x, we see that ri"''(l — x) are orthogonal polyno- 
mials with respect to the probability measure a{l — x)"~^l]o,i[(x) (ix. The 
corresponding moment sequence is 

(39) 



and the corresponding orthonormal polynomials are -yj {a + 2n)/a ri"'*(l— x). 
Therefore 

ir„(x, y) = j2 ^^ri")(l - x)T^:\l - y), (40) 

k=0 

showing that the coefficient to x'^y^ in aKn{x,y) is an integer when a G N. 
This yields 

Theorem 4.2 Let a G N. The inverse of the matrix 

1 



(41) 



a / / 



has integer entries. 

It is not difficult to prove that 



x^ 



k=0 
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and it follows that the entries of the inverse of fl41l) are given as 

n 

A;=max(i,j) 

This formula holds of course for any a > 0. 

The results of this section for a = 1,2 have been treated in the survey 
paper written in Danish. For a = 1 the formula for the elements of the 
mverse of H^n'^ was given in [B], but goes at least back to Collar [7J, while 
the formula for its determinant goes back to Hilbert in [9]. In this case the 
polynomials rn\x) are the Legendre polynomials for the interval [0,1], cf. 
[3l Section 7.7]. These polynomials have succesfuUy been used in the proof 
of the irrationality of C(3)- For a = 2 we have (a + 2k) /a = 1 + fc, so the 
coefficient to x^y^ in ( HOl) is an integer. In this case Theorem 14.21 can be 
sharpened: The inverse of the matrix (l/(^^2^"'))o integer coefficients. 
This result is also given in [14j. 

Added June 2007 A result equivalent to Theorem 12. II is given by Collar 
in [7]. Denoting by 

Mn = (pij), <i,j <n 
the matrix of coefficients of the orthonormal polynomials, i.e. 

n 

Pi{x) = ^PijX^, 2 = 0, l,...,n, 

j=0 

where Pij = for i < j, then the orthonormality can be expressed as the 
matrix equation MnHnM^ = In, hence 

= MiMn. (42) 

Collar uses f l42p to obtain formula f l5B]) and states: "Equation fH2]) . which 
provides an elegant method for the computation of the reciprocal of a moment 
matrix, is due to Dr A. C. Aitken. The author is grateful to Dr. Aitken for 
permission to describe the method and for many helpful suggestions." 

The paper by Collar is not mentioned in Choi's paper [H] and was not 
included in the list of references in the first version of this paper. 

In |2] the authors have defined a g-analogue of the Hilbert matrix for any 
complex q different from the roots of unity and have proved a g-analogue 
of §6^. When g = (1 - V5)/{1 + ^5) one can recover the results about 
the Filbert matrices and for q = —e~'^^,6 > results of Ismail ([ITj) about 
Hankel matrices of generalized Fibonacci numbers. 
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